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Observation of cooperative Mie scattering from an ultracold atomic cloud 
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Scattering of light at a distribution of scatterers is an intrinsically cooperative process, which 
means that the scattering rate and the angular distribution of the scattered light are essentially 
governed by bulk properties of the distribution, such as its size, shape, and density, although local 
disorder and density fluctuations may have an important impact on the cooperativity. Via mea- 
surements of the radiation pressure exerted by a far-detuned laser beam on a very small and dense 
cloud of ultracold atoms, we are able to identify the respective roles of superradiant acceleration 
of the scattering rate and of Mie scattering in the cooperative process. They lead respectively to 
a suppression or an enhancement of the radiation pressure. We observe a maximum in the radia- 
tion pressure as a function of the induced phase shift, marking the borderline of the validity of the 
Rayleigh-Debye-Gans approximation from a regime, where Mie scattering is more complex. Our 
observations thus help to clarify the intricate relationship between Rayleigh scattering of light at a 
coarse-grained ensemble of individual scatterers and Mie scattering at the bulk density distribution. 
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At first sight Rayleigh scattering from point-like parti- 
cles and Mie scattering from extended objects are quite 
different phenomena. For example, Rayleigh scattering 
exhibits resonances due to the internal energy structure 
of the particles, e.g. atoms, while Mie scattering shows 
resonances induced by the boundary conditions, the tar- 
get imposes to the incident light field [TJ [2] . On the other 
hand, any extended object, e.g. a dielectric sphere, is an 
assembly of microscopic scattering particles. This raises 
questions concerning the description of Mie scattering at 
objects with coarse-grained density distributions. An- 
other interesting question is, whether small ensembles of 
particles exhibit scattering features normally attributed 
to Mie scattering. 

While Mie scattering is obviously a cooperative process 
involving a macroscopic particle distribution, in quantum 
optics Rayleigh scattering is often described as a single 
particle effect. That this concept is erroneous in the case 
of ensembles of scatterers, has been recognized as early 
as 1954 by Dicke [3]. Rather, collective effects play a 
dominant role even on the level of single-photon scatter- 
ing 0H7J provided, as has been demonstrated in a recent 
experiment [8], the coarse graining of the cloud's density 
distribution is negligible. 

The scattering of light is usually studied via its angu- 
lar distribution. However if the target is an atomic cloud 
with temperature below the recoil limit, the modification 
of the atomic momentum distribution due to radiation 
pressure can conveniently be mapped by time-of-flight 
imaging techniques and may yield valuable information 
about the nature of the scattering process. Momentum 
images are routinely analyzed in experiments with ultra- 
cold or Bosc-condensed atoms, however mostly in regimes 
where local atomic disorder dominates over global Mie 
scattering. In this letter, we present an experimental ap- 
proach emphasizing the role of Mie scattering. The idea 



consists in using ultracold atomic clouds compressed to 
volumes so small, that the whole cloud acts itself as an 
inhomogeneity at which light is diffracted like at a very 
small dielectric sphere. Furthermore, instead of analyz- 
ing the whole momentum distribution, which is blurred 
by the frequent occurrence of interatomic collisions in the 
dense cloud, we concentrate on its first moment, i.e. the 
center-of-mass acceleration under the influence of radia- 
tion pressure. We discover, as a clear signatures of Mie 
scattering, an up to three-fold increase of the cooperative 
radiation pressure with respect to single-atom radiation 
pressure. Exploiting the fact that the index of refraction 
of an atomic cloud can be tuned over huge ranges by 
changing the cloud's density and volume, or by tuning 
the light frequency, we are able to address the Rayleigh- 
Debye-Gans limit of small phase-shifts for the incident 
laser beam or to enter the Mie regime of large phase 
shifts. 

Force on the center-of-mass. — We start with a brief 
reminder of the model of collective radiation pressure ex- 
erted by a low intensity incident laser beam on a cloud of 
N two-level atoms [7J |S] concentrating on cases of large 
push beam detunings, Aq 3> T, where multiple scattering 
and collective saturation can be neglected. In this limit, 
the modification of the radiation pressure F c acting on 
the center-of-mass of the ensemble of atoms with respect 
to the radiation pressure F± acting on isolated atoms can 
be expressed as [8] , 
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The averages s N = (\S N \ 2 )e, v and f N = (\S N \ 2 cos9)e tV > 
of the structure factor Sjv(k) = J2f=i e( k ~ k °) ' rj are 
taken over the total solid angle of emission of a pho- 
ton with wavevector k, at an angle 9 with the incident 
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wavevector ko [7J. They account for the role of the ge- 
ometry in the photon absorption and reemission process, 
respectively. As stated in Ref. [5] , they are well approx- 
imated by sjv — 1/N + Srx, and Jn — foot where Soo 
and foo are calculated for a smooth density distribution, 
i.e. neglecting disorder and coarse-graining in the atomic 
cloud. The averages Sqo and have been calculated 
in Ref. [7J for ellipsoidal Gaussian density distributions, 



n(r) = n Q exp[-(x 2 + y 2 )/2a 2 - 
<x r = a/k and the aspect ratio rj 



l /2(Tg], with the size 
a z /a r . 



Further approximations. — In order to simplify the in- 
terpretation of Eq. we restrict ourselves to regimes 
of negligible multiple scattering, AA^/T 2 > b n , where 
bo = SN/a 2 is the resonant optical density. This condi- 
tion implies the inequality 

4Ag/r 2 > N Soo , 



(2) 



since for large (cr S> 1) spherical clouds Soo 

(2er)~ 2 and for elongated clouds SoV 1 " 1 < Soo -1 ^ [7j- Us- 
ing the inequality |2]), we can approximate the radiation 
pressure modification ([I]) as 
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The fraction in this equation describes the role of the 
superradiant enhancement of the scattering rate T — > 
NTsn- The fact that the fraction is obtained from the 
fraction in Eq. ([I]) by substituting the factor sjv by 
signifies that, in the single scattering regime and for 
large dctunings, Ao 3> I\ disorder does not influence the 
superradiant acceleration of scattering significantly. In 
contrast disorder does shape the scattering through the 
bracket in Eq. ([3]), which represents a geometrical factor 
accounting for the angular distribution of the reemitted 
light. This bracket is responsible for Mie scattering, as 
we will see in the following. 

Fig. [ija) shows the dependence of the radiation pres- 
sure modification on Ao and N. The multiple scattering 
regime to be avoided lies in the narrow gap separating 
the regimes of red and blue detuning. 




FIG. 1: (color online) (a) Radiation pressure modification as 
a function of Ao and N. Here a = 10.4 and n = 6.3. (b) 
Dependence of Soo — foo on a and n. 



Reduction and enhancement of radiation pressure. — 
For large cloud volumes, a 3> 1, the inhomogeneity rep- 
resented by the shape of the cloud is negligible, so that 



reemission occurs predominantly into forward direction. 
Similarly, elongated clouds, r\ ^> 1, bundle the emit- 
ted light into forward direction. Both features result in 
s oo — foo — 0. Hence, the geometrical factor is dominated 
by disorder, and the radiation pressure modification is 
well approximated by just the fraction of Eq. (|3j). In 
this regime, the role of collectivity is limited to a super- 
radiant enhancement of the decay rate T. In that case, 
the radiation pressure F c /F± drops monotonically with 
increasing N. Hence, for extended clouds there is only 
cooperative reduction of the radiation pressure with re- 
spect to uncorrelated radiation pressure. Measurements 
of this effect have already been published in Refs. [5]. 

In order to get noticeable cooperative radiation pres- 
sure enhancement, two conditions must be fulfilled. 
Firstly, the fraction in Eq. ^ has to be close to unity, 
2A 3> ATsoo, and secondly the bracket must be larger 
than unity, N(soo — foo) > 1. The first condition is true 
for sufficiently large dctunings or small atom numbers. 
In this case, the radiation pressure modification is well 
approximated by just the bracket of Eq. ^ . The depen- 
dency of Soo — foo on o~ and rj is plotted in Fig. [ljb). 
The figure shows that smaller cloud sizes yield larger 
s oo~ foo- Furthermore, for a strong radiation pressure en- 
hancement, spherical clouds are more advantageous than 
elongated ones. The bracket of Eq. ^ describes the role 
of cooperative Mie scattering and its interplay with dis- 
order. In the regime where the bracket is considerably 
greater than 1, Mie scattering overwhelms the impact of 
disorder. The observation of an increase of the radiation 
pressure above the single-atom value thus constitutes a 
clear indication of Mie scattering, which has never been 
observed before. 

The relative importance of the two effects, superra- 
diance or Mie scattering, can be adjusted via judicious 
choice of the push beam detuning, A , with respect to 
NTsoo- This allows us to address both limits in the fol- 
lowing measurements. 

Experiment. — In principle, reducing the volume of a 
trapped gas to sizes not much larger than an optical 
wavelength is just a matter of cooling and compression. 
However, the compression of atomic clouds to very small 
volumes is not trivial, because the increased density leads 
to high three-body collision rates. Furthermore, below 
a critical temperature, the cloud crosses the threshold 
to Bosc-Einstein condensation and develops a repulsive 
mean- field, which resists compression [5]. 

We prepare a 87 Rb atomic cloud in the \F,mp) = 
1 2, 2) hyperfine state in a magnetic trap and evapora- 
tively cool it down to quantum degeneracy. Typically 
we have N = 10 4 ..10 6 atoms at temperatures between 
T = 100.. 1000 nK. The atom number is controlled by 
manipulating the initial atom number before applying 
the evaporation ramp. The cold atomic cloud is now 
loaded into a crossed beam dipole trap by ramping up the 
light fields within 100 ms, and then slowly ramping down 
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FIG. 2: (color online) Scheme of the experiment showing the 
dipolar trapping beams, the push beam, and the imaging 
beam. The directions of radiation pressure, a, and gravity, 
g, are also shown. 



the magnetic trap. The horizontal beam of the crossed 
dipole trap is generated by a fiber laser at A r = 1080 nm. 
The power of the horizontal beam is P r = 2 W and its 
waist w r = 25 /im. The corresponding trap depth is 
U r = ^jsr^Ih — ks x 250 fiK and the trap frequency 

is u) r = AU r /mw^. = (27r) 1900 Hz. The vertical beam is 
generated by a titanium-sapphire laser at A z = 820 nm. 
Its power is P z — 230 mW and its waist w z = 100 /im. 
The corresponding trap depth is U z ~ ks x 90 fiK and 
the trap frequency u> z — (2tt) 300 Hz. The aspect ratio 
of the trap is thus n = 6.3. 

The temperature of the cloud after evaporation de- 
pends on the atom number. However upon transfer into 
the dipole trap, the cloud is non-adiabatically heated to 
a temperature of T ~ 2.5 /tK, which is almost indepen- 
dent from the atom number. For that reason the cloud 
is thermal. From this, we estimate a cloud radius of 
a = k^/ksT/mujf, ~ 10.4. 

To perform the measurement, we now apply along 
the weak trapping (axial) direction lo z a er + -polarized 
light pulse of duration tq = 20 fxs and intensity Iq — 
95.. 730 mW/cm 2 detuned from the D2 line by A /2tt = 
0.5. .4 GHz. Immediately after the pulse (within 100 /is) 
the trap is switched off. The cloud falls in free expansion 
for T to f = 20 ms before we apply an imaging pulse in 
transversal (radial) direction. 
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Observations. — Fig. J3ja) shows a typical time-of- flight 
absorption image reflecting the momentum distribution 
of the trapped atomic cloud. With N — 5 • 10 5 atoms 
we expect high densities of no ~ 4 • 10 15 cm~ 3 and 
large collision rates of 7 co /; ~ 25000 s _1 . I.e. during 
the time tq every atom undergoes almost one collision, 
which smears out and broadens the momentum distri- 
bution. This does however not affect the center-of-mass 
momentum imparted to the atomic cloud. 

We compare the center-of-mass positions of absorption 
images of the atomic cloud with and without push beam 
sequence and determine their relative displacement Az CI7l 
[see Fig. [3Tb)]. We estimate the collective radiation pres- 
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FIG. 3: (a) Typical absorption images taken on a cloud re- 
leased from the dipole trap without (upper half) and with 
(lower half) push beam sequence. Note that about 16% of 
the atoms are lost under the action of the push beam, (b) 
Vertical integration of the images shown in (a). The natter 
curve, which corresponds to the cloud after the push beam 
sequence, is shifted towards the right by an amount Az cm . 



sure from 
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and the single-atom radiation pressure from 
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The calculated iV-dependence of the cooperative mod- 
ification of the radiation pressure force acting on the 
center-of-mass of a trapped cloud, Eq.([3]), is compared 
to measurement in Fig. [4](a) . The TV-dependence is flat 
in the disorder-dominated regime of low atom numbers. 
For larger atom numbers, the force increases first, ex- 
hibits a maximum and drops for large atom numbers. 
For increased detunings A , the maximum shifts toward 
larger atom numbers and becomes more pronounced. For 
Aq = (2tt) 0.5 GHz the experimentally accessible atom 
numbers, N = 10 4 ..10 6 , are in a regime, where radia- 
tion pressure is strongly reduced by up to a factor of 
10 [see red solid curve and symbols in Fig. |4|a)]. For 
Aq = (27r) 4 GHz we expect an increased radiation pres- 
sure, i.e. F c /Fi is larger than one, which is confirmed by 
the data [see blue solid curve and symbols in Fig. Efa)]. 

During the push beam irradiation, a fraction of typi- 
cally 20% of the atoms is lost due to hyperfine pumping 
into the F = 1 state. This is caused by a misalignment 
of the push beam with respect to the magnetic field di- 
rection resulting in a deviation of the polarization from 
perfect <r + . The resulting uncertainty in atom number is 
illustrated by the horizontal error bars in Figs. |4j Also 
shown in Figs. [4] are the curves representing the separate 
contributions of superradiance and Mie scattering. Obvi- 
ously, the data corresponding to Ao = (27r) 0.5 GHz are 
better explained by the impact of superradiance, while 
the data corresponding to Ao = (27r) 4 GHz can be in- 
terpreted as being due to Mie scattering. 

For the push beam the atomic cloud represents a dielec- 
tric medium whose mean refraction index can be approx- 
imated by n r f — 1 — <j)/<7 + ib/2a, where the off-resonant 
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FIG. 4: (color online) (a) Double-logarithmic plot of the 
measured (symbols) and calculated (lines) iV-dependence of 
radiation pressure. The red symbols and lines correspond 
to the detuning Ao = (2n) 500 MHz and the intensity 
Io = 95 mW/cm 2 . For the blue symbols and lines, Ao = 
(2tt) 4 GHz and I = 730 mW/cm 2 . The red and blue solid 
lines show calculations based on the full radiation pressure ex- 
pression Q without adjustable parameters. The red and blue 
dotted lines represent just the fraction of the expression. The 
dash-dotted lines representing just the bracket coincide for 
both values of Ao. The horizontal error bars result from un- 
certainties in the measurement of the atom number and also 
from the fact that, during the push time, to the atoms are 
pumped into a dark hyperfine state. The vertical error bars 
include an assumed uncertainty of 20% for the push beam 
intensity Jo- (b) Same radiation pressure data as (a), but 
plotted as a function of the push beam phase shift <j>. In this 
diagram the dotted lines representing the fraction coincide for 
both values of Ao. 

optical density is b — b T 2 /(4:A 2 + T 2 ) and the phase 
shift = £>A /r. Plotting the radiation pressure force as 
a function of the push beam phase shift [see Fig. Qb)], 
we infer that the data taken at large detunings are in 
the so-called Rayleigh-Debye-Gans regime of Mie scat- 
tering, characterized by small phase shifts <f> < 2ir. The 
data taken at small detunings are in the Mie regime of 
large phase shifts, <j> > 2tt. However in this regime, the 
observed rapid decrease with <j> of the radiation pressure 
is dominated by superradiant enhancement of scatter- 
ing and not by Mie scattering. This is reflected in the 
absence of a maximum in the ^-dependence of the radi- 
ation pressure. The appearance of a maximum for large 



detunings suggests a strong relationship to the first Mie 
resonance [HI H2] , which needs however to be confirmed 
by a more detailed theoretical model. The model, which 
is presently under investigation, would have to include 
the phase shift for the push beam induced by the cloud's 
refractive index. We expect that this model should also 
be able to predict the occurrence of high-order Mie reso- 
nances, which could be experimentally observed in even 
more compact atomic clouds. 

Conclusion. — We have presented measurements on ul- 
tracold atoms proving that, in the limit of the approxima- 
tion made, the expression ^ correctly describes the ra- 
diation pressure force on atomic clouds. More generally, 
it applies to extended objects that can either be ensem- 
bles of scatterers like homogeneous, ordered or disordered 
atomic clouds of arbitrary shapes and volumes, or macro- 
scopic objects like dielectric spheres. This expression 
thus represents a bridge between microscopic Rayleigh 
scattering and macroscopic Mie scattering. Neverthe- 
less, at high densities, near-field effects may come into 
play, calling for further corrections to the present simple 
model. 

Atomic clouds with small atom numbers basically rep- 
resents a randomly distributed bunch of scatterers, whose 
intrinsic disorder spoils cooperativity. The radiation 
pressure is then well described by the single-atom value. 
At large atom numbers, the atomic cloud forms a smooth 
density distribution characterized by an almost perfect 
cooperativity. In this regime, depending on the laser de- 
tuning, the cooperativity is either dominated by super- 
radiant acceleration of the decay rate, for which case we 
experimentally observed an up to 10-fold reduction of the 
radiation pressure. Or it is dominated by Mie scattering 
leading to an observed up to 3-fold enhancement of the 
radiation pressure. 
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